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one exceptional set of curves, obtained at the United States Model Basin,* 
which shows two maxima : a phenomenon which has not received explana- 
tion. It is conceivable that this may be a case in which the two maxima 
indicated in the intermediate curves of the present paper have become 
prominent through some unusual features of the model. In this con- 
nection, it must be remembered that the present calculations are based 
upon a surface pressure of specially simple type, one symmetrical round a 
point ; one could extend the calculations by integration, as in the previous 
results for deep water, so as to apply to a pressure distribution, giving 
a better analogy with ship form. It may be anticipated that the results 
would be of the same character in general, though no doubt better agree- 
ment could be obtained in certain details. 
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Introduction. — The problem of the diffraction of electromagnetic waves 
by a perfectly conducting sphere is of interest both from the point of view 
of wireless and from that of physical optics. Two cases may be considered : 
(1) when the source of the waves is a Hertzian oscillator on the surface 
of the sphere ; and (2) when the waves are plane. The formal series solution 
of both these problems has been given by several writers, including Sir 
J. J. Thomson, the late Lord Eayleigh, and Prof. H. M. Macdonald. For 
a sphere of which the radius is small compared with the wave-length the 
series converge rapidly and are suitable for computation, but for a large 
sphere the important terms are far on in the series and the latter must 
be transformed in order to get formulae which may be of use. 

For case (1) this problem has been attacked by Macdonald, Poincare, 
J. W. Nicholson, A. E. H. Love and several others, but case (2) has not 
attracted nearly so much attention. This paper, then, deals with case (2). 

* D. W. Taylor, loc. cit., p. 115. 
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The most recent paper on the subject is by Dr. Bromwich,* who gives the 
series solution in the following form : — 

At a large distance, r, from the sphere, of radius a, the scattered wave is 
given by components in the directions of r, 9, cf> respectively 

X = ca. = 0, 

3M 1 8N 



Y = Cry 



d0 sin 6 d<p 



r , 1 3M , M 

sin 9 deb d9' 



where 



M - -sin 9 cos 6 ^ 2 (-1)* -^\r |4tS IV (cos 5), 

N= ^sin^sin^~I(^l)-4 !i ^l4M^ / (cos^ (1) 

^ Ar i x ; n(n + l) ¥j n (ha) v 7 

The incident wave is of unit amplitude, is travelling along the axis 9 = 0, 
is polarised in the plane <f> = |7r, and has wave-length 27r/&. The time- 
factor e* Ac * is omitted throughout. 

S w (fe&), E w (ha) are the Kiceati-Bessel functions which have been tabulated 
for integral values of n and ha up to 10 by the British Association 
Committee, f 

Dr. Bromwich (loo, cifc) has transformed the series M and N, for ha large, 
by an application of Kelvin's " Method of stationary phase," and gets the 
approximate formulae 

Y = cos <j> (e~ ikr jhr) \ha exp (2iha cos J 9), 

Z = — sin <£ (eT ikT jhr) \ ha exp (2iha cos | 9). (2) 

These formulae, the equivalents of which were published in the earlier 
papers of Mcholson, are for points not near to the axis 9 = it behind the 
sphere; the method does not seem to be capable of dealing with points 
near to this axis, nor does it seem suitable for obtaining a second 
approximation. 

For points near to the axis 9 = 7r, Dr. Bromwich has applied another 
method, depending on the use of Green's theorem, and obtains 

Y = iY 2 cos <£ (e~ ikr /hr), Z = iZ 2 sin (e~ ikr /hr) } (3) 

where Y 2 = Z 2 = — AaJi (ha sin 0) /sin 0. 

These points are not dealt with by Mcholson. 

* Bromwich, 5 Phil. Trans., 5 A, vol. 220, p. 187 (1920). See also Nicholson, ' Proc* 
London Math. Soe.,' vol. 9, p. 67 (1910) ; vol. 11, p. 277 (1912). 
t Doodson, etc., « B. A. Eeport,' p. 13 (1914) ; p. 39 (1916). 
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With a view to verification of these formulae, numerical results* have 
been calculated directly from the original series for ha = 10. Formulas (2) 
were found to give quite good agreement for the range 6 = 0° to 90° ; for 
the remainder of the range there was marked divergence. Formulas (3) 
again give fair agreement for the range 170°-180°. 

A paper by Prof. G. N. Watsonf on the problem of case (1), mentioned 
above, suggested another method of attack. Watson transforms his series 
into another depending on the zeros of a Bessel function, and shows that 
this latter series converges very rapidly for large values of the argument, 
the first term alone giving a sufficient approximation. 

It very soon appeared, as indeed might have been anticipated from the 
form of Dr. Bromwich's approximation, that it was not possible to transform 
the series M and 1ST into others in Watson's manner, the contribution to 
the contour integrals from a large semicircle not tending to zero as the 
radius tends to infinity. 

An example in Whittaker's 'Modern Analysis'! suggested a rather 
different procedure, whereby the series M and N" are transformed into 
series depending on the zeros of Bessel functions, together with a contour 
integral along a "path of steepest descent." The contribution from the 
series is very small, and thus an approximation is given by the contour 
integral. This can be evaluated, after the manner of Debye,§ as an asymp- 
totic series ; the first term gives Bromwich's formulas (2), and a second 
approximation can be obtained. 

Comparison with the figures of Proudman shows an improved agreement 
over the range 0°~90°, but no better agreement over the remainder of the 
range. 

The method breaks down as 6 approaches ir for two reasons : 

(1) Laplace's approximation for the Legendre function P s ( — cos 6) of high 
order $ is no longer valid when 6 is so close to ir that s (ir — O)* is small. 

(2) Debye's simple approximations for Bessel functions of large argument 
are no longer valid when the order and argument are nearly the same. 

The first fact shows itself in the series which we neglect for 6 not near 
to tt; this is found to involve exponentials with negative indices Aa?*(0— 7r), 
A being a positive constant, and the terms will no longer be negligible 
when a?*(7r— f?) is small (x is written for ha). 

This difficulty can be evaded by the use of Mehler's approximation for 

* Proudman, Doodson and Kennedy, ' Phil. Trans.,' A, vol. 217, p. 279 (1917). 

t Watson, < Boy. Soc. Proc.,' A, vol. 95, p. 83 (1918). 

% 2nd edition, p. 145, Example 7. 

§ Debye, ' Math. Annalen,' vol. 67, p. 535 (1909). 
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the Legendre function in terms of a Bessel function of zero order, and we 
find that the series M and N" are now transformed into 

(a) A series depending on the roots of a Bessel function, this series being 

still negligibly small ; 

(b) Several comparatively simple contour integrals which are negligible 

in the earlier case, but important now— approximations to these can 
easily be obtained ; 

(c) A contour integral involving Bessel functions of order s and argu- 

ment x along a curve in the s-plane which passes through the point s =.x. 

The evaluation of (c) is a matter of considerable complication and labour. 
In the first place, approximations to the various Bessel functions for nearly 
equal values of order and argument are required ; some of these are given 
by Watson,* but others had to be worked out on similar lines. 

These approximations involve Bessel functions of orders + %, ± §, and the 
leading terms in (c) are found to depend upon definite numbers arising 
in the form of integrals of which the integrands are complicated expressions 
in these. To calculate these numbers, Dinnik*sf Tables of the Bessel 
functions involved are of use, but they contain several misprints, and are 
not sufficiently detailed for small values of the argument, so that a con- 
siderable amount of arithmetic has been necessary. 

The final result, for 6 near to tt, is complicated ; there is a term of 
order x 2 , which is like Dr. Bromwich's but not identical with it, and there 
are further terms in x^ z . Much better agreement is found with Proudman's 
figures for the range 170°-180°, but again the range of validity is little 
increased. For larger values of Jca, the next terms in the approximation, 
involving presumably x 2 ^ 9 would be important ; the labour necessary to 
obtain these by this method would be very great. 

The intermediate range, 90°-170°, is thus still without adequate treatment. 

(1) The series to which approximations are required are 



p—ikr oo 

M = —sin 8 cos <£-— — 2 (- 

^ kr i 


"1) 


p—ikr qc 

1ST = —sin 6 sin 6 — : — S (- 

T kr i 


-1) 


ith the notation of Watson,^ 




S» (x) = y\r n (x) 





n 



n 



2» +1 |^_) 
n(n+l)TZ n (ka) v J 



E n (x) = — i% n (x), 

where rj n (x) = (|- 7ne)*H n+ i (1) (x), g, (x) = (J wa?)*H n+i < 2 > (x), 

* Watson, < Proc. Camb. Phil. Soc.,' vol. 19, p. 96 (1916). 

f A. Dinnik, < Archiv der Math, und Phys.'(3), vol. 22, p. 226 (1914). 

% Watson, < Koy. Soc. Proc.,' A, vol. 95, p. 83 (1918). 
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and thus, since when n is a positive integer 

P n ' (-cos 6) = (-l)» +l P n ' (cos d), 
: iA. cos <£ (e~ ikr Jkr), N = ^B sin (e" ifcr /lcr) } 



we have 
where 



M 



A = S -^J^ P.' (-cos 0) sin 6 ±±M 



B»i -T^ P/ (~cos 0) sin *^, 
i n(n+l) v ' g^a?)' 

and x is written for ha, 

(2) Following the example of Prof. G-. N. Watson,* we shall consider 

Bessel functions of which the order s is a complex quantity, and we shall 

use the transformation 

s = x cosh 7 = x cosh (a -f i/3), 

whereby the s-plane, slit along the real axis except between the points ±x, 
is represented by the strip of the 7-plane between /3 = and /3 = 7r. 

Fig. 1 represents the s-plane divided by the curves into regions 



- T 




Fig. 1. 
* Watson, loc. cit. 
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(numbered 1-7) in which the Hankel-Bessel functions are given by different 

contour integrals, as shown in the Table on p. 93 of Watson's paper. 

We shall require to apply the method of " steepest descents " to integrals 

involving 

exp [2^{sinh7~- 7 cosh 7 -I- 1^ {it— 0)eosh 7}]. 

Suitable curves must pass through the point where this index is stationary, 
i.e., the point 7 = \ (iT—0)i ) and must have 

I {sinh 7— 7 cosh 7 + J ^(77-— ^) cosh 7} = const. == cos^-0. 

This gives 

cosh a sin/3— a sinh u sin /3-f- {J- (tt— #)— ft} cosh a cos ft = cos |-0. 

There is a double point at a = 0, ft = f (7r — 0), the branches there having 
tangents a = ± {/3—|- {ir—Q) }. 

Along the branch touching a = —/8 + £ (tt— 0), the real part of 

sinh 7 + { I f (?r — 0) i — 7} cosh 7 

is negative; this branch goes off to infinity at (00 , 0) and ( — 00 3 it). 

The corresponding curve, which we will call C, in the s-plane, passes 



S -plane 




y ~ plane 



JCi 




Fig. 2. 
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through the point s = x sin \Q of the real axis and goes off to infinity in 
directions parallel to the real axis but not asymptotic to it. These curves 
are shown in fig. 2. 

(3) Now consider the contour integral 

4^ sin P,- 4 ' ( - cos 0) % a ~ h 9*)'*}. T A (31) 

taken round a circuit in the s-plane consisting of any curve C symmetrical 
with regard to the origin and passing through it from the third quadrant 
to the first, together with a large semicircle below C The curve C must 
not cross the boundary between the regions 6 and 2 and between 5 and 7 ; 
it must, in fact, behave at infinity like C. 
Now in Watson's regions 5, 7, 3, 2 we have 

H.« (*) = S,« (a) -e- 2s " S s < 2 > («), S,W(a ! )- e -^S.q>-( iB ) t 

S S «(V), S,«(a;)-S,«(aO; 

respectively, where when 1 7 1 . is not small and x is large, 

S,< 2 > (#) ^e-^^v-yoosh^+iiri^^/^^^g^ (~^ 7 )}. 
&-* (0) S s < 3 > (a)' S s <*> (o^) - S.W (oj) ' 



S.W (a?) Sfi) (x) 



respectively, and 



(|7r^ S s (1) (x) I — I sinh 7 1 ~"£g p , 
(|'^ S s (3) (x) j -^ I sinh 7 1 ~%e 



— p 



where P is the real part of x (sinh 7 — 7 cosh 7). 

If s = E* 4 '*, being a positive angle not greater than ^tt, we have when E 

is very large 

sinh 7 j -— K/#, 

P — E [cos <f> (log 2K/a>- 1) - <j> sin 0]. 

Thus in 5, P is negative and of order E cos <£ log 2E/#. 

As we move round the semicircle through 7 to 3 and 2, P increases, 
taking the value — JEtt on the imaginary axis, vanishing on entering the 
fourth quadrant, and becoming positive in the remainder of 7 and in 3 
with the approximate value E cos cf> log 2E/#. In 2, P is negative like 
— E cos <]E> log 2B/#. 
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Thus 



is of order 1 on the semi-circle, except in the latter part 



&-* 0) 
of 7 and in regions 3 and 2, where it is of order e ~ 2 ' p ' . 

Also the approximate value of P s _$' ( — cos 6) sin 6 for large values of s is 
(2s/7rsm #)Ksin fs(7r— 0)— £tt}, when 6 is not near to or7r. 

Thus approximately | P s _|' ( — cos 0) sin 6 \ ^(R/2tt sin 0)*eB<*-tf)sin* Hence 
the modulus of the integrand behaves like 

(2E/7rsin #)ig-R (*■+#) sm <£ j n re gi on 5 an( j the first part of 7, 
like (2E/?r sin 0f e ~M*+e)sm4>~-2 \ p i j n tne remainder of 7 and in 3, 
and like (2U/tt sin #)^- 2 1 p I in region 2. 

Thus the integral round the large semi-circle tends to zero as its radius 
tends to infinity. 

Next consider the poles of the integrand which lie within the contour. 
These consist of the zeros of % s „±(x), and of the points s = ft+^, where n 
is zero or a positive integer. 

1 2w + l 



The residue at s = n + \{n > 0) is 

1 



. . sin0P n '(-cos0)£^. 



The residue at s = A is 



2 



^ sin ,^ ir cot A 0. 
2r- 



:7r 



As regards the zeros of J s _| (a?), Watson* has shown that these are simple 
and not on the real axis, and in fact lie near to the dotted curve proceeding 
from the point s = x into the fourth quadrant. 

The residue at a zero s = v is 



^sinflP,_/(-co S fl) „^iM 



v'-i 



[8/a S T s -f(«)] s =,« 2w +i' 

Thus we get, after some re-arrangement, 
B = —ismxe %x c,ot\d 

+ 2iri 2 ^-r sin ?„_.'( -cos 6) ^.^r^fl, rr— r 

+ 1 7=1 sin * • p -»' (- 008 *> • tjlr ?=r+I & 

(4) Next consider the contour integral 



(3-2) 






VH sin (9 P s / ( -cos 6) I s - 1 J4 277 i *, 



(41) 



taken along the curve C. 



In region 5, 






6 -2E7rsm^. j n re gi on 2, 
* Watson, £oc. ciV., § 5. 



ifr-l 0*0 



1, and thus 
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the modulus of the integrand at each end of the curve behaves like 
(K/27rsin0)^~ R(w+a)sin ^ and the integral converges. 

Now we have H-,< 1 >(oj) = e S7d M. s M(x), H- s < 2 >(¥) = *r 5 «H s < 2 > (x), and thus 

' * 9 ) { -z—. — - = I* *) ( — T~- — t> % ' e -> tnis expression is an even function 

of s. So also is P 8 -b (—-cos 0). 

Hence the integrand is an odd function of s and the integral vanishes. 
Adding to the equation (3*2), we get 

•B = 



ism%e ix Goi^0 



+ 2m 2 vj(v*—\) sin ?„-/"(— cos0) 



^-j(a?) 



(^•*+i)p/&(;. i (»)] t » l 






I O 1 ^* 



sin . "P s -h (—cos 0) 



'2±s~k (a?) 



1 , Vs-ki®) 



w§» 



(4-2) 
The integral in equation (4*2) can be written 

'2i|r,_j(a:) 1 



J. 

1 









sinflP^i' (—cos 0) 



sin0P 5 _| / (--eos0) 



- Hr) 



e 



■2sTri 



"%_*(*) , 1 



+ U 



&. 



ds (4-3) 



(4-4> 



The form (4*3) is suitable on the upper half of the curve, each part con- 
verging at this end ; the form (4*4) is similarly suitable on the lower halt 

Equation (4*2) has been proved for not equal to or tt; it can easily 
be seen that it holds for = tt, but it does not hold for = 0, as the 
Legendre functions of non-integral order become infinite in this case. 

(5) As regards the integral I& we note that sin 0P s _ i / ( — cos 0) can be 
expressed as the sum of two terms 

where Ci and G2 can be -expressed either as definite integrals or as asymptotic 
series in s~ l not involving exponentials.* 
Now consider the integral 

2s 



T ' 



S* — f fs-fW 



the circuit consisting of C and the large semi-circle below C. 

On the part of the semi-circle below the real axis, the integral tends to zero 
on account of the factor g-«(«— *)*; on the part above the real axis (in 

region 2) it tends to zero on account of the factor r*-fw 

* Watson, * Canib. Phil Trans,/ vol, 22, pp. 291, 295 (1918). 
VOL. 0. — A. 2 L 
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Thus 



v 






js-i 



\ x )\s-v 



(5-1) 



The remaining part of I2 is 

s f 2ty s ^(x) 






C/S 2_l 



' grgfr-fll-frr]^ _ S111 iZizl ( ~ C ° S ^^ ^ (5*2) 







— 257Tt 



+ 1 



To get this in a form involving an integral along a " curve of steepest 
descent," we must consider the appropriate approximations for ty s ^(x) and 
&-$(#) at various points of the curves C and C. These curves pass through 
Watson's regions 5, 1 and 2, in all of which £*-h(w) is represented 
by (|w)*S s (2) (x); but in which t; 5 _x(^) is represented respectively by 

(iwx)*{8 t to(m)-e- M 8 9 ®(x)}, (Jww)*S,W(a;) and (Jt^)*^) (a)- S, «(*;)}. 
Hence, as the integrand for which is a suitable curve must contain the 
factor B s ^(x)/S s (2) (x), which of course is discontinuous at the points Q, S, 
where C crosses the boundary curves between the various regions, the integral 
involved must be (see fig. 1) 



•Q 



P s 2 -i 



e [s(«-9)-lK}iQ 1 



L?~*(x) 



+ #* 



'2srri 



+ 



qs 2 -J &-*0*0 



s s 2 -~i 



s 






~T* X 



e?s. 



The remaining part of I 2 " is, on reduction, 



-f 



p <S — -g 



e -iH*+o)+WGids- 









c* 



~%L 



pi . ^ p 



# 



2s?n 



-f" J- 



e 



— 2srri 



+ 1 



efo. 



This last integral can easily be seen to vanish over the large semi-circle 
below C, and thus the contribution of this last term is 



00 






sin 0T?n ( — cos #)— J cot J0. 



1 n(n+l) 

The series can easily be evaluated; its sum is —J tan §#. This contribution 
is therefore equal to — cosec 6. 
Thus 

cosec 0- 



t " 
±2 






»T 









f £ s 



S/ 1 ) (x) 



tS,^(a?) 



^N^-H^Ci^ 



2 ^ S -^-r s [-"('-^)+Hi c 2WT57 ^^ — . (5*3) 



This is for not equal to or ir. 
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(6) Collecting the results so far obtained we get 
B = -J tanJ0-Je*»cot $e + 2iriS v/(v 2 -i) . Vv ^(x)f[dfds^(x)] s= , v . 



C^- _.._. , . C 2 






e 



• 2«7rt 



_j_ j^ 



•Q 



-i T «[.(M-wc 1 (fe + 
s « 2 -i 



P * — 4 

s S f < 1 >(a!) 



e -[.(»+i)+WiC l( ls 

1 



° ^ ¥4 g[s(7r-fl)-iir]* Q ^ /A.I \ 



Precisely similar work applies to the series A ; the only difference being 
that we have r) s _^' (x) and &-%' (x) in place of rj s ^(x) and £,_$(#), and 
consequently that the residue at s = |- which arises in the equation 
corresponding to (3*2) gives a term cos % . e ix . cot -J# . 

Thus 

A= -Jtan|0 + ^ 2 -cot§0 + 27n^ 



r-/,[v(7r-6l)-|7r]i fi—[v{ir — G) — \n\i . " 
° p v pi 



2v7Tt 



+ 1 



e 



■2v7rt 



+ 1 



_ o<e J. •*• 

P * X 



•T 



S 5 4 



e [«(ir-fl)-l-tC 1 & + 



c s 2 — \ d/dx{xiS s ( 2 Ux)} 1 

(6-2) 

(7) The next step is to investigate the co-ordinates of the points Q and S, 
with a view to determining the magnitudes of the terms in equations (6*1) 
and (6-2). 

By actual plotting of the curves it is found that the imaginary co-ordinate 
of Q lies between 06 # and 1*5 &, while the imaginary co-ordinate of S lies 
between 0'6x and 0. 

Hence when x is very large, it is seen that the integral from P to Q 
involves an exponential, the real part of whose index, (jr + 0) times the 
imaginary co-ordinate of Q, is large and negative. This integral is thus very 
small. The same is true of the integral from S to T, except when 6 is near 

tO 7T. 

In this latter case, S is given very approximately by 

The index of the exponential involved is thus — ^~ . (?r — 0f } and thus the 

integral is no longer very small when (V — 6) is of the order x~*. 

(8) We have next to consider the terms in A and B (equations 6*1, 6*2) 
arising from the zeros of £,_$' (x) and &_.$(#) respectively. 

For &_j' (x) the work is given by Watson,* whose results are as follow. 

* Watson, ' Koy. Soc. Proc.,' A, vol. 95, pp. 90-97 (1918). 

2 L 2 
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The large zeros of $,_$' (x) are given approximately by the equation 

#(sinh 7-- Ycoshry) — \iri = — (m-\-^)m, 
where m is an integer large compared with x, and the corresponding value 
(v m ) of v which makes % v _% (x) have a simple zero is 

Vm ~~ (m -*- i) 7r 2 / 2 (log w) 2 — ( m + 1) irij log m. 
The corresponding value of 7 is 

Jm ^ — log 'W' + I TT^. 

These zeros lie in Watson's region 7, in which 

^(a?) = (i7^x)^S s (V(x)-e~ 2s « i S s W(x)}/(l^e^ s «% 

Using the known approximations for S s (1 >(#) and S s ®(x) when |y| is not 
small, we find that at the zero z^ OT 

Vv~h( x )l$fi s &-*'(#)],=„ — l/2y m ~ - 1 / 2 lo g m - 
Hence since the imaginary part of z^ m is negative, the term in A arising 
from a large zero of &_$' (#) is approximately 



7T% 



r e ~(m + l)ir(w + 0)f\ogmQ^^ e - (m+l)ir{Tr~ff)l\ogm Q^ 



vm log m 
which is very small as long as 9 is not nearly equal to ir. 

For the small zeros of &-*'(#)> Watson* obtains the approximate values 

where for the first three 

/> = 0-8083, 2-577, 3-88. 

Using these results, we arrive at the conclusion that the first terms in the 
series in equation (6*2) contain exponentials with indices 

—~h\/3p • (w + 0)#* and ~~ 1 4/3/0 (tt — 6) x% 
and axe thus again small, for large values of x, except when is nearly equal 

tO 7T. 

The work for the corresponding terms in B (equation 6*1) is very similar 
and gives similar results. As these terms are negligible, it is unnecessary to 
put down the details of the calculations. 

(9) We have shown that for 6 not near to tt (more accurately, as long as 
(tt — 6)x* is large), we have approximately, neglecting terms containing 
exponentials with large negative indices, 

(9-1) 



B = -itanii?-ie 2ix cot^ + 



c s 2 ~~i S s < 2 >(^) 



S ^ (1) (4 g[.(»-«-Jir]< 0l ^ (9-2) 



* Watson, foe. c^., p. 97. 
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It remains to approximate to the integrals along the curve C. We must; 
first write down the asymptotic expression for Ci. For the corresponding 
part of P 5 _i(--cos#) we have, |sj being large, 

5 [ S (^e)-^]t(27r5sin 0)-*/l + -L cot 01 



Differentiating, we get 



3i 



Oi^ -i(«/2vsin^)*(l-^cotd). 



(9*3) 



The asymptotic expansions for S s (1 > (x) and S s (2) (x) are given by Debye* ; he, 
however, uses ir where we have 7. 

So long as 1 7 1 is not small, which is true along the curve C when is not 
near to 7r, we have 

Sd) 
(x) ~~ exp {x (sinh 7 — 7 cosh 7) — J nri} { \ irx sin ( — ^7) } ~* 

1 + ( o ""• 04 C0th2 7 ) x S * nh 7 

S(2) 
(^)-^exp { — #(sinh 7— 7C0sh7)4-j7r^} {\irx sin (—^7) }~s 

1— (- — — coth 2 7J /^sinh7 

After some reduction we get for the integral of equation (9*1) ? neglecting 
terms of relative order x~~ 2 , 

(2tt sin 0)~% s~~i exp {2x (sinh 7—7 cosh y) + ix cosh y(ir— 0) — \nri] 



and for equation (9*2), 



1-5 icotd 4- 1/4 + 7/12 co thy 
8 x cosh 7 x sinh 7 



ds (9*4) 



(27rsin#) M s * exp {2# (sinh 7—7 cosh 7) + ^ cosh 7 (tt— 0) — \m] 



c 



1 _3 frcotfl 1/4-5/12 coth 2 7 ' 
8 a? cosh 7 x sinh 7 



ds. (9-5) 



The curve C is, as we anticipated, a "path of steepest descent "for our 
integrals, the " saddle-point " being at 7= \i {nr—0). 
To evaluate them, we put 

2sinri7~ 27 cosh y + i(7r — 0) cosh 7 == 2icosJ0 — r 2 . 

Along C, t 2 is real and positive, and the integrals involve e~ xr<i multiplied 



* Debye, ' Math. Annalen,' vol. 67, p. 535 (1909) ; ' Munchen Sitzungsber.,' Abh. 5 
(1910). 
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by a factor which may be expanded in positive integral powers of t. The 
individual terms, of the type 



J —oo 



exp( — XT 2 )r n dr, 



are known integrals, and thus the evaluation may be effected to any term in 
the asymptotic series, though it becomes extremely laborious even to include 
terms of relative order x~ 2 (for which, of course, additional terms are 
required in the approximations, both for the Legendre and the Bessel 
functions). 

The final results obtained are 






• J tan ^0 -f | e 2ix cot -| + 1 cosec J exp (2ix cos J 0) 

{1 + \i sin 2 \ 9 sec 3 \ 0x~ l 
- i (5 tan 2 J + 7 tan 4 \ 0) sec 2 J 0aT 2 } 



r<&T 2 }, (9'6) 

B 



= — ■ £ tan |- — | e 2i> cot \ — J cosec | exp (2£e cos \ 0) 

{1 — \i sin 2 J sec 3 \0x~ x 
+ i (5 tan 2 10 + 7 tan 4 J 0) sec 2 J tor 8 } 



Now Yi + iY 2 = i 



Zi + iZg = i 



~dA B_; 

.30 sin 0. 

dB_ A " 
.90 sin 0_ 



(9*7) 



On differentiating, it is found that the terms arising from the first two 
terms in equations (9*6), (9*7) vanish. 

The values of the quantities that we are investigating are thus 
approximately, for not near to tt, 

Yi-HYs = \x exp (2ix cos J 0) {1 + \ix~ x sec 3 \ 0— \ x"~ 2 sin 2 J sec 6 § 0}, 

Zx + ^Z 2 = — J a5 exp (2ias cos J 0) {1 + J i»~ 1 cos sec 8 J + f- ^?~ 2 sin 2 -J sec 6 J 0}. 

The following Table shows the values of Y ls Y 3 , Z h Z 2 calculated for 
ha = 10 from the second approximation above. These are denoted by the 
letter A; and for comparison the corresponding values obtained by Proudman, 
Doodson, and Kennedy are also given under the letter B. 
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ha = 10. 



e. 



o 

10 
20 
30 
45 
60 
70 
80 
90 
100 
110 



Y x . 



B. 



1 
2 
3 
4 
4 

•3 
•4 

4 
3 



Y 9 . 



A. 



•76 


2 02 


•16 


2'2L 


•12 


2-yo 


•34 


4*26 


•78 


4-98 


•60 


0-20 


•62 


~3'27 


•84 


-4-41 


•73 


-1-40 


•53 


4-51 


•47 


4-22 



4-67 
4-52 
3 -92 
2-51 
1-52 
•4-98 
•3-48 

1-37 

5-00 
2-33 

•3-84 



B. 



4-38 

4-57 

4-16 

2-30 

-1-35 

- 5 -09 

-3-83 

1-96 

5*09 

0-97 

-3-27 



Zi. 



■1 

2 

•3 

■4 

•4 

■0 

3 

4 



•4 

■1 



•76 
•16 
•13 
•36 
•76 
•40 
•81 
•67 
•02 
•84 
•89 



B. 



-2-02 
-2-32 
-3-16 
-4-31 

-4-79 

-0-40 

3-85 

4-72 

0-14 

-4*92 

-2-45 



Z 2 . 

A. B. 



4*67 

•4-52 



3 

■2 
1 

4 



•91 

•47 
•61 
•99 



3-25 



1 
5 



•80 
•00 



1-28 
4-65 



-4-38 

-4-36 

-3-93 

-2-51 

1-64 

4-98 

3-34 

- 1 -75 

-5-19 

-1-52 

4-81 



(10) As a preliminary to considering values of near to 7r, let us take the 
case in which = tt. 

In this case we must differentiate our series before evaluating them. We 
have from §1 : — 

Y = -cos $ l - '-j-t ( - 1)» 4^_ 

§4 £ {sm 0P,' (cos*)} -liP,' 
.K act E» 

! ~ ar i( iy 2w+1 

^r i ' t&(t& + 1) 

S n 3 



(cos 0) 



Z 



sin (f> 



|-^{sinW(cos0)}-|<P,/ 



Now 



™ {sin 6 P n ' (cos 0) } = ^ (ti + 1) P« (cos (9) 



■|=r P.' (cos 0) 

COS0P n '(cO3#). 



1(101) 



Also, when = tt, (-l)*p n ( C os 0) = 1; (-1) W P„' (cos (9) = -— J^ + l). 
Hence, for = 77% 



Y = —cos </> 







—ikr 



—ikr 



J(A + B), Z=-sin^i_-i(A + B), (10'2> 



&r 



where A and B denote the series 



To the approximate evaluation of these series we now proceed. 
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In Watson's notation 



1 bn\ x ) 



CO 



B = iS(2»+l)&^. 



By work which is exactly parallel to that of paragraphs 3, 4, 



JHL — — 



sm#e' 



t# 



2ttS 



VVv-i^) 



+ ^ 



■icosaje^-- 27r 2 



(^«+i)[a/3»&_ i ( a3 )] S!=! 



+ ~T 



0' Us-fO) ^ + 1J 






V.-*(aj) 



+ 



1 



dte, 



<&?. 



(10-4) 



The series in (10*4) gives a negligible contribution when a? is large, each term 
involving an exponential with large negative index. 

As regards the contour integral, we may note that in this case, the curve G 
as defined in §2 degenerates into the boundary between regions 5 and 7, the 
boundary between regions 1 and 3, and the real axis from x to + co . 

The expression which we shall have to integrate over C is, for A, 

The remainder, by a process similar to that of § 5, is seen to be 



Q 



se 



' 28nt ds + 



s 



/»rr 



Jp^+l 



els 



s e~ 2s7d +l 



els + 



-s 



o e 2s ^ + 1 



ds 



ise 



~2sni 



2 



7T 






1 



47T 2 



6" 



■2sTrt 



■2 



*s 



JQ Jo^^+l 



efe + I s^s 
o 



= a negligible quantity + — + hx 2 , 
since S is the point x, and 



*t 



■00 



oe 



~2stri 



+ 1 



(is 



1 



1+0*^ 



1 



1 1 

l OS 



We have, finally, to approximate to 



47r 2 1 *" 2 2 ^3 2 



48* 



c S,< 3 >(0) 



- &. 



Along the boundary between regions 5 and 7, |y| is not small, and we have 
the approximation 

/ p 4 =• exp {2x sinhy—ycoshy) — |-7ri}. 

Since sinh 7 —7 cosh 7 is real and negative on this boundary, the contribu- 
tion from this part of C involves an exponential with a large negative index 
and may therefore be neglected. 
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In region 1, approximations valid right up to the point x are :* 
S s (1) (#) = — -f i tanh 7 exp {x (sinh 7—7 cosh 7) + 77} 

S 5 (2) (as) = f i tanh 7 exp { — (# sinh 7— 7 cosh 7)— ??} 

{^^'I-i/a^ + e^Ii/B^)}, 

where 77 = %x sinh 3 7 sech 2 7 , and the phase of r\ is to vanish when 7 is real 
and positive. I is the modified Bessel function. 

Thus the contribution from this part may be taken as 

— s exp {2a? (sinh 7—7 cosh 7) + 2<n\ ... . T ' / 7 — 1/<WtT / \ **• 

The principal part of this integral, of course, arises for 1 7 1 small. We 
now express everything in terms of 77. 

After some reduction, the term of the integral involving the highest 
power of x is found to be 

_ f 3-1/3^-1/3 ^ 6 -I-l/3fa) + ^^Il/3(^) , 

Now the initial phase of 7 (at x) is §7r, and hence the initial phase of rj is 
2 77, and thus, writing rj=:te 27r \ t being a real positive quantity, the approxi- 
mation to this part of the integral is 



3-l/3^4/3 e -l/37« 



W /-l/3 I-1/3 (Q — I l/3(Q j, 

i- 1/3 (0+^ 1/3 ^ii/3(0 



In an exactly similar manner, the approximation to the part arising from 
region 2 is found to be 

3-1/3 x „z e - ljan f" r i /3 L^h^i^L (U. 

JO 1-1/3 (O + ^^Il/sCO 



Adding, the final value of 



S,W(a?) . 

^ ri /o\ / \ '^* 1° 

,. C b s W(^) 



j 



Q-l/3 r 4/3 ,-l/3ir* f°° /-1/3 { I ~l/3(0"-Il/3(0} {2L-i /3 (Q+ Il/3 ft) } ,, 

Jo I~l/3 2 (0+I-l/3(OIl/3(0-V(0 

= 3-173^/3,5-1/3^ p ? gaj> 

In the same way the approximation to the corresponding integral in B, 
s 57^ — c V9 a 7iw \ f ^ 1S found to be 

_ Q-l/3 r .4/3 ,-l/3*i f 00 /"^ { J-2/3 (Q — 12/3 (0 } {1-2/3 (0 ± 2Ia/g (Q } ,, 

Jo I-2/8 3 (0 + 1-2/8 (0l2/8(0 + VW 

= -S-VS^/Sg-l^Q^ sajt 

* These approximations are not given explicitly by Watson, but can be obtained by 
application of the methods of his papers (loc. cit., § 8, and f Proc. Camb. Phil. Soc,,' 
vol. 19, p. 96 (1916) )• 
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The two integrals P and Q are definite numbers which can be calculated, 

to a sufficient approximation, by actually plotting the functions. For this 

purpose the tables of the modified Bessel functions of orders + ^-, + f , as 

given by Dinnik,* are of use, but additional figures, for small values of the 

argument, are also needed. These have been calculated directly from the 

series. For extremely small values of t, and for large values of t, the Bessel 

functions may be replaced by appropriate approximations, and thus it is only 

for intermediate values of t, say from t = 0*02 to t = 2*2, that the curves 

need be plotted. 

The final results are 

P = 1-39, Q = 1-24, 

which are probably accurate to the second place of decimals. 

(11) Collecting together our results, we have for the leading terms 

A = \ ix 2 + 1-39 x 3-V3 e 1 ^ x^ 
B = \ ix 2 - 1 -24 x 3 - V3 e i/6« cc m 

It is useful at this stage to compare the results obtained with those 
arrived at by direct arithmetical evaluation. In the paper of Messrs. 
Proudman, Doodson, and Kennedy, to which reference has already been made, 
the various series involved are calculated for x = ha = 10. It is possible to 
pick out what corresponds to the above series. In fact, the series A above is 
the same as the C + iD of equations (6), 2 of the paper quoted, and B is 
C' + iD'. 

From Tables V-VIII, we conclude that for x = 10, 

A = 19-58 + 259-64:, B = - 21-00 + i43'48. 

The approximations above give 

A = 17*98 + i60'38, B = - 16-04 + i40*76, 

which is as close agreement, 10 not being very large compared with unity, as 
can be expected. 

When, however, we come to take the sum of A and B, the terms in x^ z , 
being very nearly equal and opposite, practically cancel, and the terms 
arising from the next approximation, involving presumably x 2 ^ } will become 
important. We find 

| (A + B) = |-^ 2 +0-08 + 3- 1 /^ 1 ^^/3 ) 

* A. Dinnik, < Archiv der Math. u. Physik/ (3), vol. 22, p. 226 (1914). Some mis- 
prints may be noted. The correct values of I-i/sO^)* 1-2/3 0*0 are as f°how : — 

X. -■-— 1/3W* X I_2/3W- 

2-2 2*5626 2*8 3*7595 

3*0 4*7754 8*0 415'01 

6*0 66-55 
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and, in the notation of the first paragraph, 

Y 1 = Z 1 = -0-045 (fo) 4 / 3 , 
Y 2 = Z 2 = -£(fa*)*-0-026 (kafi*. 
For ka = 10, these give 

Yx = Z x = -0-97 ; Y 2 = Z 2 = -50-57, 
whereas the calculated values are 

Ti = Z x = +0-71 ; Y 2 = Z 2 = -51-56. 

It would thus be desirable to carry the approximations to our contour 
integrals at least one stage further, but the labour involved seems to be 
prohibitive. 

(12) There remains to discuss the case of nearly equal to tt. It is more 
convenient to differentiate the series before transforming them. 

n—ikr 

We get Y = —cos <j> -~ — [2 2 — cos 0X 4 — 2 3 ], 

rCT 

p—ikr 

Z = —sin 6 — — [2i — cos S3— 2A 

kr 



GO 



where 2 X = 2 (- 1)* (2» + 1) P M (cos 0) |4t^. 

1 ili„ (/ot) 

S 8 = I ( - 1)- (2« + 1) P» (cos 0) |4|g, 

2 ^!(- 1 ) s ^TI) P "' (co ^ ) W 

The work of transformation is exactly parallel to that of the previous 
cases. 

For 2i we have finally to approximate to 

S P s _ i (-cos^)|^i|& 

taken over the same curve C as in the case of = tt. 
The remainder is 



fS 

4- 5P s _j(— cos 0)ds. 



Of these thefirst may be neglected. Also S is the point s = x. 
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Since is near to 7r, we may replace the Legendre function in the third 
integral by J (2s cos |0), and the integral becomes 

f K 
s J (2s cos |- 0) ds = §# sec -| Ji (2as cos -J #), 

Jo 

which tends to \x 2 as # tends to 77*. 

Also the important part of the second integral is from | s | small, and so, 
swinging the contour round into the imaginary axis, we get 

3 £J (2£cos|<9) 



2 



di. 





l + e^- 



This can be evaluated as a series in cos 2 J#, the coefficients involving 
Bernoulli's numbers. The first terms are 

For the integral over C the process is practically as in paragraph 10, save 
that we have an additional factor J (2s cos §• 0). 

Thus the result will be 

0-964 Jo (2x cos \ 8) e 1 ^ 1 x^. 

Altogether the leading terms of % ly S2 are 

2i = Jiajsec JflJi(a»cosi5) + 0"964J (2iBoosJ^)« 1 / <w a^, 

v 2 - 1 i x gee \ 0J X (2x cos \ 0) -0*86 J (2a? cos \ 0) e 1 ^ 1 x^ 3 . 
In the series S3 the extra integrals are 

1 f K s 

P,_j'(-cos 0)^~ -jds-h J o ^rrj P,_j' (-cos 0)<fe, 



JoS 2 -i 


P 


/(- 


which the first 


gives 


1 
48 sin 2 £ 





1 
160 



cos 2 J #4- terms of higher degree in cos \0. 

The second integral can easily be shown to be approximately 

cosec 2 { 1 — Jo (2^ cos 1 0) } . 



The integral 



T P *-A~ C0S #)q%7^ <fc gives, on putting for 



S 2_i -«-* v wu "'s,< a >(a-) 
P 5 _^ ( — cos#) the approximate value J 5 sec -J0 cosec 2 J#Ji (2s cos J 5), the 
expression ^- cosec 2 -J (9 sec J Ji (2x cos -| 0) x 0*964 e' 1 ^ x^ s . 

Thus for the leading terms 

2 3 = — ^cosee 2 #{l — J o (2a?cos|0)} 

- 0*482 cosec 2 J sec J Ji (2a; cos | (9) e 1 /^ ^ 
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and in a similar way 

2 4 = — ^cosec 2 0{1 — J (2^cos -|#)} 

+ 0*43 cosec 2 -| sec J J x (2x cos J 0) e l l M x 1 ^ 

To compare with the figures for x = 10, we note that 

2i = C + «D, 2 a = C' + iD', 2 3 = A + iB, 2 4 = A' + £B', 
in the notation of Proudman's paper. 

The calculated values of the series give, for = 170°, 

Si = 5*65 + 37'47i 2 2 = -512 + 28-84i, 2 3 = ~7'08-24'20i 

2 4 = 6-25-18-04i, 
while the above approximations give 

2i = 6-7 + 37'2i 2 a = - 5-9 + 29*9^ 2 3 = - 6*03 - 24*3i 

2 4 = 5-38-17-7i. 

The following Table shows the values of the series Si, 2a, 23, 2 4 , calculated 
from the above approximations for ha = 10 in the range of values of from 
170° to 180°. The final figures for Y b Y», Z h Z 2 (denoted by A), given 
therefrom by the formulae 

Yi + iY 2 = 23 — cos #2 4 — 2 a , 
Zi + iZ 2 = 2 4 + cos 23 — 2i, 

are also shown, and for comparison the values obtained of these by Proudman 
(denoted by B). 



e. 


-5j. 


IV 


-*3 


• ' 


s 4 . 


170 
172 
174 
176 
178 
180 


6-7 + 37 -2i 
10 -24 + 44*69$ 
12 -50 + 50 *69$ 
15 -86 + 56 "05$ 
17 '45 + 59 -35i 
17 -98 + 60 -38i 


-5*9 + 29-9* 
-9 -13 + 33 '51$ 
-11 -94 + 36'58i 
- 14 -15 + 38 '72i 
-15*56 + 40*29* 
-16-04 + 40-76^ 


-6-03- 

-7-oi- 

-7-84- 
-8*46- 
-8*87- 

2- 


-24 -3i 

-26*30$ 
-28*27$ 
-28-7K 
-29-76* 

Si 


5*38-17*7$ 
6-25-18-64i 
6-99-19-39$ 
7-55-19-89$ . 
7-91-20'OOi 


6, 


Yi. 


Yo 


z, 


* 


rr 
Zi 3 . - 


-Ol» jD. 


A. B. 


A. 




A. B. 


170 
172 
174 
176 
178 
180 


-5 3 -7*12 
-4-07 -4-76 
-2-87 -2*33 
-1*84 -0-23 
-1 *21 -60 
-0*97 0*71 


-36*5 -35-28 
-41-35 -4L*0 
-45-57 -45*6 
-47*59 -49*3 
-49-99 -51-0 
-50-57 -51-56 


4-5 
2-95 
2-29 
0*13 
-0-68 
-0 97 


6-59 

4-7 

3 

1-6 

0-8 

0-71 


-30-9 -31-67 
-39-29 -38-0 ! 
-41-96 -43*9 1 
-47 30 -48-4 
-49*68 -51-0 
-50-57 -51-56 j 

i 



